An experimental route to spatiotemporal chaos in an extended ID oscillators array 
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We report experimental evidence of the route to spatiotemporal chaos in a large ID-array of 
hotspots in a thermoconvective system. Increasing the driving force, a stationary cellular pattern 
becomes unstable towards a mixed pattern of irregular clusters which consist of time-dependent 
localized patterns of variable spatiotemporal coherence. These irregular clusters coexist with the 
basic cellular pattern. The Fourier spectra corresponding to this synchronization transition reveals 
the weak coupling of a resonant triad. This pattern saturates with the formation of a unique do- 
main of great spatiotemporal coherence. As we further increase the driving force, a supercritical 
bifurcation to a spatiotemporal beating regime takes place. The new pattern is characterized by 
the presence of two stationary clusters with a characteristic zig-zag geometry. The Fourier analysis 
reveals a stronger coupling and enables to find out that this beating phenomena is produced by the 
splitting of the fundamental spatiotemporal frequencies in a narrow band. Both secondary instabil- 
ities are phase-like synchronization transitions with global and absolute character. Far beyond this 
threshold, a new instability takes place when the system is not able to sustain the spatial frequency 
splitting, although the temporal beating remains inside these domains. These experimental results 
may support the understanding of other systems in nature undergoing similar clustering processes. 

PACS numbers: 05.45.Xt, 05.45.Jn, 47.52. +j, 47.20.Ky 
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I. INTRODUCTION 



Understanding complex systems in physics, chemistry, 
biology and neuroscience is a topic of theoretical and ex- 
perimental relevance, particularly when we are dealing 
with the route to spatiotemporal chaos in spatially ex- 
tended systems. In this route, the existence of various 
subsystems or domains with spatiotemporal synchroniza- 
tion, coherent patterns, allows to handle a description in 
terms of clusters. Accordingly, we find cluster dynam- 
ics in networks of coupled oscillators such as: reactive- 
diffusion systems in chemistry [l|, Q , multicellular tissues 
in neurobiology [!, |j[ and coupled arrays of oscillators 
in physics 0, ||| . Although the scientific community has 
kept a great attention to complex dynamics during the 
last decades, a universal route to chaos in extended sys- 
tems has not been already described, not even under- 
stood, in contrast with the theoretical and experimental 
effort devoted to confined systems. Our aim is to char- 
acterize experimentally the dynamics of a lD-spatially 
extended system in the weak turbulence regime where 
a spatiotemporal splitting is produced inside synchro- 
nized clusters. Theoretically, Hohenberg and Shraiman 
approached the understanding of spatiotemporal chaos 
with the generalized Kuramoto-Sivashinsky equation Q 
and also there are several recent numerical works behind 
this purpose in the field of complex networks [1, Qj3] ■ 

The experimental results reported here are enclosed in 
the field of extended, dissipative and complex systems 
involving many degrees of freedom. Our experiment is 
a quasi-lD thermoconvective system, an array of non- 



linear coupled oscillators consisting of ascending convec- 
tive cells (hotspots) exhibiting a spatiotemporal beating 
regime of the type of a zig-zag pattern (ZZ). It is charac- 
terized by the presence of stationary clusters, developed 
from a previous mixed pattern of irregular clusters. If 
we consider that strong nonlinear couplings between am- 
plitude and phase might be neglected, the ID-array is 
displaying "phase locking phenomena" [ill ]. We try to 
approach a phase description to follow the dynamics far 
from the threshold of the primary bifurcation. A par- 
ticular point of view is also tackled on the phase syn- 
chronization theory of N-coupled oscillators which was 
firstly introduced by Winfree [l2l HU and later devel- 
oped by Kuramoto [lj]. Modelling extended ID-system 
of N-coupled oscillators undergoing global phase insta- 
bilities under certain constraints, such as weak coupling 
between oscillators subjected to the same driving force, 
has allowed to build rich phase spaces [15[ . 



ID-experiments are particularly interesting for the di- 
versity of patterns they exhibit, for example in: ther- 
moconvection [H, fijl 111 [H HO, El 1H, electro- 
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convection [2J, |25[, the printer's instability [26|, l27j . 
directional solidification [281 ] . directional viscous finger- 
ing [29j], circular liquid column arrays [3(| [M, 
the Taylor-Dean system [33| . Theoretical approxima- 
tions to ID-oscillatory instabilities have already been 
successfully modelled for supercritical bifurcations with 
the Kuramoto-Shivashinsky phase equation (KS) [341 ] . 
and for subcritical bifurcations with a Ginzburg-Landau 
equation adding a quintic stabilizing term [35l [3a ] . Other 
ID-localized patterns such as solitons, and spatiotempo- 
ral intermittency, have been modelled using a modified 
Swift-Hohenberg equation [3^ and the KS [38j |. In a 
previous work [391 ] we had characterized quantitatively 
primary and secondary instabilities for the same exper- 
imental cell. In this research we reported the existence 
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of a mixed state for thin layers from a stationary cellu- 
lar pattern (ST). This mixed pattern (ST/ ALT) consists 
of localized oscillatory domains in an alternating pat- 
tern (ALT) coexisting with the basic cellular pattern, ST, 
which remains always stable. Depending on the values of 
the control parameter, the new unstable pattern shows 
a different pinning to the underlying stationary cellular 
pattern ST, and thus, there might be to two different 
kinds of fronts: fluctuating and stationary. 

Mixed states in weak transitions to chaos have 
been usually identified as spatiotemporal intermittency 
regimes i.e. the Rayleigh-Benard convection in an annu- 
lar gap [13] and the Faraday experiment [4l|. A similar 
behavior is found in experiments with open shear flows 
where the basic state is always stable, such as plane Cou- 
ette [HJ], circular Couette and Poiseuille [43[. These in- 
stabilities have in common a "subcritical branch" sent to 
infinity [3]. In this sense, bistability in these systems 
is understood as the coexistence between "patches" dis- 
playing the new pattern and the original one. 

Experimentally, certain extended systems in closed 
flows show a higher complexity stage before achiev- 
ing spatiotemporal chaos or defect mediated turbulence: 
regimes with temporal or spatiotemporal frequency split- 
ting. For example, the dynamics of a f D-front in electro- 
convection with liquid crystals @, [H[ shows time biperi- 
odicity, and in the Taylor-Dean system (33j spatiotempo- 
ral biperiodicity appears. In 2D systems, spatiotemporal 
chaos via a zig-zag instability has been exp erimentally 
observed in liquid crystals [HI, [46|, H3, 51] and it has 
been successfully modelled in the frame of the Ginzburg- 
Landau equation for a weak nonlinear regime [4^ | and for 
a selective forcing [5fjj . 

According to the results reported here, we have no 
knowledge of any previous experimental work of a zig-zag 
regime (ZZ pattern) in a lD-convective system caused by 
a spatiotemporal beating, and developed from an irreg- 
ular clustering process. So, this report is the first ex- 
perimental evidence in hydrodynamics of an array of os- 
cillators displaying clusters, irregular and stationary, as 
the result of an increasing coupling interaction between 
hotspots. 

We report a cascade of secondary bifurcations in a 
quasi- ID convective system, a rectangular layer of fluid 
opened to the atmosphere under a ID-heating at the cen- 
ter of the bottom plate. For a fluid with Prandtl num- 
ber (= viscous diffusivity/ thermal diffusivity) Pr 75, 
the competing effects between buoyancy and capillarity 
account for a Benard-Marangoni convection mechanism. 
For a fixed depth [d = 3 mm), as we increase the con- 
trol parameter (the vertical temperature difference AT V ), 
the system undergoes a subcritical instability towards 
a mixed pattern of irregular clusters from a stationary 
cellular pattern. Further on, a spatiotemporal beating 
regime settles down, this consists of one or two large sized 
clusters with stationary fronts. The quantitative analy- 
sis by complex demodulation techniques shows a contin- 
uous evolution of the selected order parameter, the am- 
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FIG. 1: Top view (a) and cross section (b) of the experimental 
cell. Both sketches do not keep the same scaling. 



plitude of the critical modes, from which the spatiotem- 
poral beating regime is found to be supercritical. Both 
oscillatory instabilities to the mixed pattern (ST/ ALT) 
and to the beating regime (ST/ZZ), have absolute and 
global nature. We will focus on the behavior of the criti- 
cal modes nearby the thresholds of these transitions and 
describe the subsequent lost of spatial frequency splitting 
of the oscillatory modes (ST/DW) in a new transition. 



II. EXPERIMENTAL SETUP 

A fluid layer of depth d is placed in a narrow rectangu- 
lar cell L x x L y (L x = 470 mm, L y = 60 mm) [Fig. [TJa)]. 
This layer lies over a mirror [Fig. [T{b)] in order to im- 
plement the shadowgraphy technique. At the center and 
along a line in the x direction, the layer is heated at 
Th by means of a thermoregulated rail of 1 mm (thick- 
ness). This rail is placed beneath the mirror and it is 
laterally isolated providing a smoothed gaussian temper- 
ature profile in the y direction . The lateral walls are 
two aluminum blocks (coolers) whose temperatures are 
kept constant at T c = 20.0 ± 0.1 °C by means of a sec- 
ondary water circulation. A more detailed description of 
the experimental setup is found in [39| . 

The fluid used is silicone oil with a viscosity of 5 
cSt. This fluid is transparent to visible light with a 
slow evaporation rate, therefore in the range of temper- 
atures explored we can work with the free surface of the 
layer opened to the atmosphere whose temperature is 
controlled at T a = 22.0±0.I°C. The control parameter 
for a fixed depth is the vertical temperature difference, 
AT V — Th — T a . The relevant physical properties of the 
5 cSt silicone oil are: the thermal conductivity A — 0.117 
Wm^K -1 , the thermal diffusivity n =6.68 I0~ 8 m 2 s _1 , 
and the surface tension a = 19.7 mN m . 

The results reported here are measured for a fixed 
depth d — 3 mm. The most suitable adimensional num- 
ber for the Benard-Marangoni convection mechanism is 
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FIG. 2: Instantaneous shadowgraphies of the patterns: (a) ST 
for the stationary mode M 3 (k s , 0); (b) ST/ ALT for the oscilla- 
tory domains with counter propagating modes M v ±(k v , ±w v ) 
which are resonant with M s ; (c) ST/ZZ for the spatiotempo- 
ral beating regime with new modes M z ±(k z ,u) z ) very close to 
the previous M v ±. Below we show a zoomed sketch of the 
oscillators involved and their characteristic wavelengths (A s , 
A-y, \ z = A^ <5A). 



the dynamic Bond number Bop = R/M (=Rayleigh 
number/ Marangoni number), for which at AT V = 30 
K we get Bob ~ 1-04. In consequence, thermobuoyancy 
effects (R) are balanced by thermocapillary effects (M). 
There are two basic characteristic time scales: the vis- 
cosity time scale t v « 2 s and the thermal diffusivity 
time scale t k 135 s, hence the convective time scales 
are supposed to be subordinated to the diffusive temper- 
ature scales. 

The shadowgraphic flow- visualization technique allows 
us to observe thermal gradients in the bulk of the fluid 
layer when an incoming parallel light is sent through the 
convection pattern. Once it is reflected back at the mir- 
ror surface, the output beam is projected into a screen. 
The space-time varying modulation of the light beam on 
the screen reveals the characteristic spatiotemporal peri- 
odicity of each pattern, in Fig. [2] we show some images. 
The darker aligned spots in the center correspond to the 
uprising thermal plumes (or hotspots) at the heating line. 
The screen image is recorded with a CCD video camera 
with a resolution of 570 x 485 pixels and constant focal 
length, aperture and gain. The optical devices allow us 
to visualize on the screen a centered area of 310 mm x 
60 mm from the total surface area (450 mm x 60 mm) 
[Fig. Ha)]. 

From the image on the screen (Fig. \2§ the acquisition 
image system records along a line two kinds of spatiotem- 
poral diagrams: one over the bright modulation parallel 
to the aligned hotspots at a fixed contrast position y Q , 
S Vo (x,t) [Fig. [3ja,d)], and another one perpendicular to 
it at a given position x , S Xo (y,t) [Fig. Ufa)]. 

The following results and theoretical discussions con- 



cern the study of a lD-convective system whose geomet- 
ric aspect ratios, at d = 3 mm, are T x = L x /d = 150 
and T y — L y /d — 20, hence the ratio jr- = 7.5 allows to 
classify the system as weakly confined in x. 

Measurement process 

Spatiotemporal diagrams S Vo (x,t) are processed via 
the bidimensional Fourier transform and complex de- 
modulation techniques in order to determine the mod- 
ulus of the amplitudes \Aj\, the wave numbers kj and 
the frequencies u>j of each fundamental Fourier j-mode 
M(kj,ujj) [see Fig. |3jc,f)]. So, each oscillating pat- 
tern will be defined by the following field: ^(x,t) = 
J2j Aj exp [i (kjX + uijt)] + c.c. Also, demodulation tech- 
niques allow us to characterize the dynamics of localized 
patterns and fronts by defining the following parameters 
on S Vo (x,t): 

(a) The signal S Vo (x, t) is processed by dismissing the 
normalized amplitude above a critical value /i = 
1/e [see Fig. [3Jb,e)]. This method is used to de- 
termine the average size of localized patterns, L c , 
and the nature of their boundaries: stationary and 
fluctuating fronts. Thus, the parameter L c is a 
coherence length. A similar parameter is defined 
in other systems like in the Kuppers-Lorz insta- 
bility (sij . the plane Couette flow [52j and the 
Belousov-Zhabotinsky reaction [53|. 

(b) The velocity of propagation of boundaries, v p , is 
determined by measuring the velocity of the front 
defined by the new traveling mode with /i = 1/e. 
This method is not suitable for fluctuating bound- 
aries. 

(c) The invasion rate, o~i nv , is the ratio between the 
total surface occupied by the domains belonging to 
a certain mode with respect to the total extension 
of Sy o (x,t). This parameter is measured for the 
new unstable mode or modes invading the whole 
system. 

(d) The attenuation length, £, is measured for coex- 
isting patterns, and it is the length for which the 
new unstable pattern penetrates further beyond the 
coexisting domain [33, [H| . We choose the unsta- 
ble mode which represents the new pattern, and 
we measure the length for which its amplitude de- 
cays according to A ~ e^ x ^ (from critical phe- 
nomena (5f| ), where £ is a correlation length of the 
system. This parameter might be understood as a 
"spatial memory" inside the neighboring domain. 
This measurement is not suitable for fluctuating 
boundaries. 

From the unidimensional Fourier transform of the sig- 
nal S Xo (y,t), qualitative and quantitative information 
can be extracted about the dynamics of a particular 
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FIG. 3: Spatiotemporal diagrams S Va {x,t) in the ST/ALT (a) and the ST/ZZ (d) regimes; (b) and (e) are the corresponding 
filtered signals by complex demodulation techniques with fi — 1/e; and (c) and (f) are the corresponding Fourier spectra where 
the fundamental modes have been defined. 



hotspot. We have checked that outside the acquisition 
position selected to record S Vo (x,t), the temporal fre- 
quencies of a given oscillator are the same than the ones 
provided by S Vo {x,t). The sequences of measurements 
are taken for ascending steps of 1 and 0.5 K, and corre- 
spond to the asymptotic states. 



III. RESULTS AND DISCUSSION 

A. A general approach to the experimental results 

In the bulk of the fluid layer, for every control param- 
eter, a primary convection (PC) appears. It consists on 
two counterpropagative rolls which ascend at the heat- 
ing line and descend at the lateral cooled walls. As soon 
as AT V ^ 0, a primary bifurcation to the stationary ST 
pattern takes place (see the stability diagram in Fig. O . 
In the shadowgraphy image of ST [Fig. UK a)] we observe 
this multicellular pattern, with wavelength A s ~ 6 mm 
[A s w 2d, see Fig. Ufa)]. The physical aspect ratio de- 
fined as Ta s = L x /\ s « 80 is large enough and thus, 
under this condition the system is extended. The ST 
pattern is represented by the stationary mode M s (k Sl 0). 
At d = 3 mm this pattern never looses stability against 
the emerging unstable traveling modes. This is the cause 
for which we always find localized domains of the new 
pattern coexisting with the ST pattern. As we increase 



the control parameter AT V , the transition to weak turbu- 
lence starts. For the different regimes, results on k and 
uj (in Fig. O and on modules \ A\ (in Fig. [7]) are obtained 
from the demodulation analysis of the spatiotemporal di- 
agrams Sy a (x,t). The subsequent secondary instabilities 
give rise to the following asymptotic patterns summa- 
rized in table [J (see also the stability diagram in Fig. [5]): 

(i) A mixed ST/ALT pattern. At AT vcl = 21 K, 
the system bifurcates towards a regime of fluctu- 
ating domains, these are irregular clusters of vari- 
able size in the ALT pattern that usually collapse 
[Fig- E^b)]. This ALT domains show a period 
doubling of the basic pattern (k v = k s /2), with 
two counterpropagative modes M v ±(k v ,±uj v ) su- 
perimposed to the stationary mode M s [k s , 0) . Ac- 
cording to the weakly nonlinear coupling revealed 
from the Fourier spectra [i.e. Fig. [31(c)], inside the 
fluctuating fronts the two traveling modes locally 
destabilize the basic ST pattern, so the resonant 
triad [M s (k s ,0), M v ±(k v , ±uj v )] is triggered. Ac- 
cording to results in Fig. Efb) traveling modes have 
an average period of the order of 15 s. 

(ii) A spatiotemporal beating regime or ST/ZZ pat- 
tern. At AT VC 2 = 31 K, and for a range of 5 K, 
one or two large domains in ZZ, with an average 
width of L c w 80 mm coexist with the ST pat- 
tern [see the coherent ZZ domain in Fig. [2Kd,e)]. 
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The ZZ domains are stationary clusters resulting 
from the splitting of the spatiotemporal frequen- 
cies of the fundamental modes. The spatial beat- 
ing phenomenon comes from the existence of two 
close wavcnumbers (fc z , k v ) of the traveling modes 
[Fig. HJd)], meanwhile their characteristic zig-zag 
geometry is due to the temporal beats [two close 
temporal frequencies (u> v , lo z ) extracted identically 
from diagrams S Xo (y 7 t) in Fig. S^b) and diagrams 
S Vo (x,t) in Fig. 0Jc)]. From the bidimcnsional 
Fourier spectrum [Fig.^f)] we confirm the increas- 
ing number of harmonic modes taking part in a 
stronger nonlinear dynamics. So, this beating phe- 
nomenon corresponds to the following fundamen- 
tal modes: the stationary mode M s (k s , 0), the new 
emerging traveling modes M z ±(k z , ±w z ), and the 
previous existing ones M v ±(k v , ±lo v ). From the re- 
sults on Fig.[B{a,b), we obtain a narrow spatiotem- 
poral frequency band (5k, 6u>), where 5k — \k z — k v \ 
and 5u> — \uj v — u) z \ are the envelope periodicities 
of the spatial and temporal beats, respectively. 

(iii) A temporal beating regime or ST/DW pattern. 
From AT vc3 = 34.5 K [Fig. EKb)] onward, the 
spatial frequency splitting disappears. So, in this 
regime the competition between the split traveling 
modes is not able to sustain spatial beats, only tem- 
poral beats remain by means of a new mode which 
becomes unstable in a unique and larger localized 
domain. This domain is composed of the follow- 
ing fundamental modes: M s (k s ,0), M v ±(k v ,±u) v ), 
M w ±(k v ,±uj z ). 

In Fig. [21 we show the characteristic distribution of 
hotspots for these regimes: ST, ST/ALT and ST/ZZ. 
Discontinuities in the dynamics of hotspots denote the 
presence of ID-fronts. We may think of hotspots as a 
sort of "coarse-grained units" linked to the hydrodynam- 
ical variables which are able to oscillate like limit-cycle 
oscillators. Each characteristic pattern depends on the 
couplings between oscillators, future work is devoted to 
determine different types of couplings, this topic is dis- 
cussed in the last section. A review of patterns for thicker 
layers is found in (39j . 

We have found no hysteresis for neither of these insta- 
bilities, i.e. there is no different behavior in the descend- 
ing sequence compared with the ascending one presented 
here. For higher depths, classical subcritical bifurcations 
have been found 1391 1, From the diag rams Sy o (x : t^, in 
Fig. E£a,b) we show the evolution of the modulus of the 
amplitude of each fundamental mode belonging to the 
same quadrant of the symmetric Fourier spectrum: M s , 
Aft,_, Ajf z _. In order to avoid large values of \A S \ due 
to the resonance of the stationary mode at the bound- 
aries [3{|, data in Fig. \T[b) is determined inside the 
boundaries confining ZZ and DW domains. 

When we disturb the surface locally and close to the 
threshold of the bifurcations to ST/ALT (at AT vcl ) and 
ST/ZZ (at AT VC 2), the corresponding oscillatory domains 
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FIG. 4: In the ST/ZZ regime, at AT V = 33.5 K : (a) Spa- 
tiotemporal diagram S Xo (y,t) for an oscillator belonging to 
the ZZ domain; (b) Fourier spectrum in uj from the signal 
S Xo (y, t) [in (a)] at the acquisition position to record S Vo (x, t). 
From the signal S Vo (x,t): (c) The Fourier spectrum in u> [the 
scaling is different from (b)]; (d) The Fourier spectrum in k. 
The components of the critical mode M z + (k z , w 2 ) and the pre- 
vious existing one M„+(fc„,w„) have been indicated in their 
respective graphics. 
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FIG. 5: Stability diagram. Continuous lines bound regions 
with the same asymptotic dynamics. Discontinuous arrows 
correspond to the critical points, the continuous upward ar- 
row corresponds to the ascending sequences of measurements. 
Stationary patterns are: PC (primary convection) and ST 
(cellular pattern). Oscillatory patterns are: ST/ ALT (mixed 
pattern of irregular clusters), ST/ZZ (spatiotemporal beat- 
ing regime of stationary clusters), ST/DW (temporal beat- 
ing regime) and TW (traveling waves). ST/ALT+TW corre- 
sponds to the coexisting ST/ ALT and TW patterns 



become globally unstable in transient regimes. Besides, 
fronts show no transversal drifting, so the group velocity 
is null and both instabilities have an absolute character. 

For consistency with the following sections, we define 
two reduced control parameters e\ — AT V / AT vc i — 1 and 
£2 = AT tI /AT! UC 2 — 1 for the two subsequent instabilities 
towards the mixed ST/ ALT pattern and towards the spa- 
tiotemporal beating regime ST/ZZ, respectively. 



B. Irregular clustering dynamics in the mixed 
pattern 

In Fig. [5] we show a sequence of filtered spatiotem- 
poral diagrams in the mixed ST/ ALT regime where the 
surrounding ST pattern has been filtered. The system 
is bistable in this regime because localized domains in 
ALT, with well defined and fluctuating boundaries, coex- 
ist with the basic cellular pattern ST. Fluctuating bound- 
aries are irregular fronts which bound these localized 
domains, hotspots belonging to these domains are syn- 
chronized in lo v and therefore the images in Fig. [5Ja- 
d) remind us the nature of irregular clusters (eventually 
they collapse) like the ones appearing in the Belousov- 
Zhabotinsky reaction-diffusion systems [ij. Above the 
threshold [e\ = 0) clusters spread over the whole dia- 
gram S Vo (x, t) [i.e. Fig. HJc)] until they stop collapsing 
at AT V « 26 K (ei = 0.24), and saturate into a unique 
domain of great spatiotemporal coherence [i.e. Fig.[8{e)]. 
This domain will remain in the following instability to 
spatiotemporal beats [i.e. Fig. [S[f)] where two new un- 
stable modes will emerge. 



The basic cellular pattern ST does not ever become 
unstable in the wide range of the explored control pa- 
rameters (AT V = [21,43] K). According to previous ex- 
perimental work (see references in pij). the bifurcation 
diagram would correspond to that of a global subcritical 
bifurcation where the threshold branch is unstable for 
a diverging control parameter e% (see sketch in Fig. [SJ. 
If we take into account thermocapillary effects compet- 
ing against thermogravitatory effects, from the dynamic 
Bond number Bojj oc d , we infer that by decreasing 
d, thermocapillary effects become stronger in the mixed 
ST/ ALT regime than for higher values of d, where this 
dynamical regime coexists with a traveling wave pat- 
tern [3^ (see the stability diagram Fig. [5]). Consequently, 
we may consider that at d = 3 mm subcriticality is sent 
to infinity, and bistability is a natural consequence in the 
following instabilities where critical modes become un- 
stable inside these confined domains. However, the sub- 
sequent instability might not necessarily be subcritical as 
it will be shown in the next section. 

In this transition, the amplitudes of the traveling 
modes are constant (see Fig. |H] where \A v+ \ ss l-A^-l) 
in the range < e x < 0.35 (e x = AT v /AT vcl - I). The 
frequency u> v increases monotonically with AT V and has 
an average value of 0.42 s _1 (< ujt v >= 0.75) [Fig.[B£b)]. 
In Fig. UOf a) we show the invasion rate of clusters ai nv 
vs. s\. Following the sequence of pictures in Fig. [5] we 
infer that, from an initial domain in ALT, as we increase 
the control parameter (ei), Oinv increases until it spreads 
inhomogeneously (pi nv ~ 0.4) all over the spatiotemporal 
diagram S Vo (x,t). This process corresponds to the left 
branch in Fig. lIOf a) with fluctuating fronts. Similarly, in 
other transitions to turbulence in open flows we may say 
that, at this stage, the system is being "contaminated" 
(i.e. plane Couette flow [42]). Once irregular clusters are 
equally distributed over the spatiotemporal diagram, the 
average size of these irregular clusters is L c « 30 mm. 
The spatial and temporal correlation lengths of the sta- 
tionary mode decay in this regime. This is the cause for 
considering ST/ ALT as a spatiotemporal chaos regime. 

From this point onward, di nv decreases until it achieves 
a practically constant value of approximately 0.14. This 
is due to the fact that we are settling the system close 
and below the following instability with the formation 
of a unique ALT domain of great spatiotemporal coher- 
ence and with stationary fronts (v p ~ 0.01 mm/s). This 
behavior corresponds to the right branch in Fig. If Of a). 
We have checked the same behavior below the following 
bifurcation with a smaller step. Regarding these results 
there is a growth of the "phase coherence" which is sup- 
posed to fit a similar role as the invasion rate when it 
increases in the left branch. In this left branch, as we 
increase the control parameter more oscillators are syn- 
chronized nearby uj v according to an irregular clustering 
dynamics. 

The theoretical approach to this clustering process is 
close to the diffusion-induced inhomogeneity proposed by 
H. Daido and K. Nakanishi [56| except for the global cou- 
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TABLE I: Regimes in the route to spatiotemporal chaos for increasing AT V 



Stationary pattern (ST) extended pattern 1 mode: M 3 (k a ,0) 

Mixed pattern (ST/ ALT) localized domains 3 modes: M s (k s , 0), M v ±(k v , ±uj v ) 

Spatiotemporal beats (ST/ZZ) localized domains 5 modes: M a {k s , 0), M v ±(k v , ±u v ), M z ±(k z , ±uj z ) 

Temporal beats (ST/DW) localized domains 5 modes: M 3 (k 3 ,0), M v ±(k v , ±w„), M w ±(k v , ±w z ) 
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FIG. 6: Evolution of the wavenumbers (a) and the frequencies (b), respectively scaled with d = 3 mm and t u = d 2 jv = 1.8 
s. Vertical dashed lines correspond to threshold values (AT VC 2 = 31 K, AT vc z = 34.5 K) which bound the ST/ZZ regime and 
continuous lines in the ST/ZZ regime are a guide to the eye. 



pling, if we suppose that hotspots are being locally cou- 
pled by diffusion and convection. The growth of phase 
coherence is supposed to trigger the clustering process 
until the "phase synchronization" spreads over the whole 
array of oscillators, when there are two stationary clus- 
ters (ST/ZZ regime). 



C. Supercritical bifurcation to the stationary 
clusters regime 

For even higher temperatures the system becomes un- 
stable towards the beating regime. In this ST/ZZ pat- 
tern, studied along two ascending sequences by steps of 
0.5 and 1 K, two domains saturate to a fixed width value 
of L c w 80 mm. The ALT pattern inside has become 
unstable because of the emerging unstable modes. These 
are two new traveling modes, M z ±, which are very close 
to the previous existing ones, M v ±. From Fig. [S]Ja,b) we 
extract, in the ST/ZZ regime, the envelope frequencies: 
5k w 0.04 mm -1 and 5lo w 0.10 s _1 . These results are in 
agreement with the spatiotemporal periodicities of the 
demodulated signals shown in Fig. [TTT a-d). the spatial 
periodicity of the beats is of 157 mm (« 2L C ) and the 
time periodicity is of 124 s (« 62 s of maximum ampli- 
tude, brightness periodicity of zig-zags). The beating dy- 
namics of the ST/ZZ pattern depends on the competition 
between the unstable modes M z ± and M v ± , their ampli- 



tudes and the spatiotemporal frequency band (Sk,Suj). 
The invasion rate is approximately constant Ci nv — 0.18 
[Fig. HM. 

Above the threshold (£2 = 0) this regime is strongly 
nonlinear and beyond this instability, the pattern looses 
the spatial frequency splitting. Above AT vc s « 34.5 K 
(s2 = 0.11), the temporal frequency splitting still remains 
in the new ST/DW pattern [Fig.[7[a,b)]. In this ST/DW 
pattern, a unique and wider domain in DW exists whose 
ID-fronts, as we increase AT V , convert into a chain of dis- 
locations [i.e. Fig. [T27a.b)]. Besides the stationary mode 
M s drifts as it has already been observed in other sys- 
tems for high control parameter values (57| . Inside this 
DW domain we observe an increasing number of tempo- 
ral beats (their period diminishes) the further the system 
is from the threshold, in agreement with the results in 
Fig. Ellb) where the distance 8u> becomes larger. Further 
above £2 > 0.29 a different dynamics due to the increas- 
ing number of dislocations is responsible for a variable 
Ui„„ [Fig. fTDTb)]. In Fig. [7Ja) we notice that the evolu- 
tion of the rate of the amplitudes \A Z _ |/|^4„- 1 increases 
continuously from 2 up to 7 until the necessary condition 
(in k) for the existence of spatial beats does not accom- 
plish at the ST/DW regime. 

The continuous growth of the modulus of the ampli- 
tude of the critical mode M z _ has been checked analyz- 
ing both the whole diagram S Vo (x, t) and the ZZ domain 
[Fig. E(a)]. This continuous evolution of the modulus of 



8 



(a) 







□ 












ST/ZZ 


ST/DW 


▼ 






□ □ D □ m B 
i 



(b) 



9 □ □ D □ 



0.06 0.12 0.18 0.24 0.3 0.36 



1-0rr 
0.9- 
0.8- 
0.7- 
0.6- 
0.5- 
0.4- 
0.3- 
0.2 
0.1- 
04: 



II 

ST/ZZ i ST/DW 







0.06 0.12 0.18 0.24 0.3 0.36 



FIG. 7: (a) Evolution of the amplitudes of the fundamental modes M z -, M v - inside the domains defined by ZZ and DW; (b) 
Evolution of the amplitude of the stationary mode M 3 inside the same domains in order to filter the resonant contribution of 
the mode M a at the boundaries [the scaling in amplitudes is different from (a)]. Dashed lines correspond to threshold values 
and continuous lines are a guide to the eye. In (a) and (b) £2 = AT V / AT VC 2 — 1. 
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FIG. 8: Sequence of filtered diagrams S Vo (x,i) with /z = 1/e, the reduced control parameter is defined as ei — AT V /AT VC \ — I 
with ATud = 21 K. Diagrams (a)-(e) correspond to the mixed ST/ALT pattern, diagram (f) corresponds to the beating regime 
ST/ZZ at the threshold (e 2 = 0). 



the critical amplitudes \A Z ±\ agrees with a supercritical 
bifurcation. On the other hand, from results in Fig. [6] it 
follows that the phase velocity — uj/k along the whole 
cascade grows linearly with the vertical temperature gra- 
dient. 

The supercritical transition to ST/ZZ is the beginning 
of a phase synchronization transition nearby the station- 
ary cluster phase determined by the critical frequency lu v . 
This phase transition has to be completed far from the 
threshold £2 = [58|. According to subcriticality (the 



subcritical ST branch is imposing bistability) , from the 
analysis of the correlation length as an attenuation length 
of the new pattern (ST/ZZ) inside the original one (ST), 
we get: < £ >ss 25 mm. This subcritical length cannot 
be determined with fluctuating boundaries, nevertheless 
for the mixed ST/ ALT pattern the role of <Ji nv has con- 
tributed to the interpretation of a sort of diverging cor- 
relation length that follows the inhomogeneous growth 
of the phase synchronization over the whole array of os- 
cillators [Fig. HPT a)]. Regarding the supercriticallity of 
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FIG. 9: Evolution of the amplitude of the fundamental 
traveling mode M„_ along an ascending sequence of ei = 
AT V I AT vc i — 1 inside the bifurcated domains. Vertical dashed 
lines separate different regimes and the continuous line corre- 
sponds to the averaged amplitudes. The inserted sketch cor- 
responds to the diagram of bifurcation to the ST/ ALT regime 
from the stationary ST regime. 



the critical modes in the ST/ZZ regime M z ±(k z , ±u> 2 ), 
the correlation length £ at the threshold is supposed to 
diverge, so its comprehension goes beyond the point of 
view of critical phenomena regarding the attenuation of 
critical amplitudes. 



IV. FURTHER REMARKS AND CONCLUSIONS 

For a thermoconvective fluid layer, we have reported 
a route to spatiotemporal chaos for a system of high- 
dimensionality with N w 80 coupled oscillators. Syn- 
chronized interactions between hotspots define domains 
with the same oscillatory pattern and with specific ID- 
fronts. These ID-fronts can be of the fluctuating type, 
for irregular clusters (with weak coupling between os- 
cillators) , which take part in a process of inhomogeneous 
growth of the phase coherence. For higher control param- 
eter values a new type of clusters emerge with stronger 
coupling between oscillators: ID-fronts are stationary in 
this phase synchronization transition to ST/ZZ. Bifur- 
cations to the mixed ST/ ALT (irregular clusters) and 
ST/ZZ (stationary clusters) patterns are global instabil- 
ities where synchronized clusters have different coherent 
widths expressed in terms of L c . We have shown that 
the size of the clusters at the beating regime is given by 
the wavelength of the spatial phase envelope and comes 
from a spatial frequency splitting. 

In order to follow the whole sequence of bifurcations 
presented here, if the primary instability from PC to- 
wards the stationary ST pattern is produced at AT V = 
0.1 K, we may define the reduced control parameter 
e = AT.„/0.1K —1. The bifurcation diagram in Fig.[T37a') 



shows the cascade of bifurcations in the transition to 
spatiotemporal chaos for this particular extended sys- 
tem. Using the amplitude of the critical modes, as the 
most convenient order parameter, we show how the sys- 
tem undergoes a supercritical bifurcation to the beating 
regime ST/ZZ. There is a subcritical width in the ST/ZZ 
pattern determined by L c which is compatible with the 
existence of a supercritical bifurcation of the unstable 
traveling modes M z ±. This instability appears from a 
mixed ST/ ALT pattern of localized clusters which are 
inhomogeneously spread. Oscillators belonging to this 
spatiotemporal chaos regime are synchronized in u> v , but 
often these irregular clusters collapse, similarly to an ag- 
ing clustering dynamics for weak coupling. The growth 
of phase coherence accounts for the clustering process un- 
til the ensemble of N-coupled oscillators becomes phase 
synchronized in lo z far from the threshold of this beat- 
ing regime ST/ZZ (AT v2 ) [Hj]. Beyond this ST/ZZ pat- 
tern the system keeps only temporal beats, ST/DW. We 
have sketched the evolution of the clustering process in 
Fig. [T37 b) where the order parameter is the module of the 
amplitude of the critical traveling modes, \A\. The dy- 
namics in the ST/ZZ regime depends on the competition 
between the split traveling modes. 

The existence of localized domains in the beating 
regimes (ST/ZZ and ST/DW) approaches the statisti- 
cal model for chaos of Hohenberg and Shraiman Q in 
extended ID-systems, when the scales of energy supply 
and dissipation arc of the same order: the depth of the 
fluid layer d (being d <C L x ). Under these conditions, at 
d = 3 mm our system exhibits the existence of two spa- 
tiotemporal coherent domains of width L c < L x . L c > £, 
where £ has been defined as a correlation length of the 
system and measured as an attenuation length of the 
new ST/ZZ pattern inside the surrounding pattern ST 
(< £ >~ 30 mm and ^ < Ij). In this context, the do- 
mains of width L c are correlated surfaces in space and 
time. Nevertheless, if we keep in mind the whole dy- 
namics (from the mixed ST/ ALT pattern to the beating 
regime) and look for a physical explanation to these kind 
of transition phenomena, we are pushed to understand 
synchronization transitions in a convective system. For 
a given ensemble of N-limit cycle oscillators distributed 
on an ID-array of extension L Xl depending on the type 
of diffusive coupling, the interaction range between oscil- 
lators, Xj can be: (i) local (between neighboring oscilla- 
tors, x — L X /(N — 1)) leading to spatiotemporal chaos; 
(ii) global (all-to-all, \ — Ac) leading to collective syn- 
chronization phenomena; and (iii) nonlocal (further than 
the nearest neighbors, L X /(N — 1) < \ < Ac) leading 
to the existence of coherent patterns. Cooperative phe- 
nomena are present in both secondary bifurcations, the 
interaction between hotspots corresponds to a diffusive 
coupling between nearest neighbors which, as far as we 
increase the control parameter, becomes stronger due to 
nonlinearities. A phase dynamics may be suitable for 
the weak nonlinear regime at the mixed ST/ ALT pat- 
tern (this work is in progress). Moreover, a phase de- 
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FIG. 10: (a) Invasion rate of the traveling mode M v - close and above the bifurcation to the mixed ST/ALT pattern for two 
different sequences by steps AT = 1 and 0.5 K. The reduced control parameter is e\ = AT V / AT vc i — 1. Continuous lines are a 
guide to the eye; (b) Invasion rate determined by the modes M z _ and M v - representing the beating regime with \jl = 1/e and 
along an ascending sequence. The reduced control parameter is £2 = AT V / AT VC 2 — 1. Vertical dashed lines separate different 
regimes at threshold values: AT vc i = 21 K, AT VC 2 = 31 K, AT vc3 = 34.5 K. 




FIG. 11: A sequence of demodulated diagrams S Vo (x,t) selecting the pair M„_, M z _ for increasing AT V in the ST/ZZ regime: 
(a) 31 K (e 2 = 0), (b) 31.5 K (e 2 = 0.02), (c) 32 K (e 2 = 0.03), (d) 33.5 K (e 2 = 0.08). Dark and bright regions correspond to 
the minimum and maximum values of the amplitude. 



scription can be extended to the beating regime (results 
on nonlocal coupling for this regime will be reported else- 
where (58J), although its strongly nonlinear character, ex- 
cept for the defect mediated turbulence displayed in the 
ST/DW regime. Regarding the Kuramoto model [l4| 
for a weakly nonlinear system of interacting i-oscillators 
(for i = 1, . . . , TV) with their respective raw i-phases, fa, 
which include the contributions from the fundamental 
and the harmonic modes, the most suitable model might 
approach the following phase equation for a local (non- 
local) coupling: 



N 

fa = F t (uj t ) + Y / Hij (*) ' (fa - fa ) (1) 
3=1 

Where Fi is a function of the fundamental frequen- 
cies of the i-oscillators: Fi(ui v ) in the mixed ST/ALT 
regime, and Fi(u) v ,u z ) in the beating regimes. Hij(t) ac- 
counts for the time-dependent adjacency matrix whose 
ij-elements are nonzero, Hij(t) = 1, for each connected 
pair of ij-oscillators, while Hij(t) = otherwise. In our 
system, this connection defines the topology of the net- 
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FIG. 12: (a) Spatiotemporal diagram S Vo (x,t) at AT V — 42 K (e2 = 0.32); (b) the corresponding diagram for the phase gradient 
of the mode M a , fronts of dislocations are observed. 
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FIG. 13: (a) Bifurcation diagram \A\ vs. e, where e = AT„/0.1K —1 is the reduced control parameter defined from the 
threshold of the instability to ST while e(AT„ c i), e(AT„ C 2) and e(AT„ C 3) are the reduced control parameter at the thresholds 
of the subsequent instabilities; (b) Sketch of the dynamics of the localized domains for the same increasing control parameter 
e, the order parameter is \A\ for the corresponding critical modes. 



work which happens to be linked to the dynamics (i.e. 
beats [581 ] . collapsing and spreading phenomena). In fact, 
Hij{t) is providing the range of interaction, meanwhile 
Tjj represents the magnitude of the coupling strength be- 
tween connected oscillators. Near threshold, usually 
takes the easiest description (501: ~ sin(0j — <j>j), but 
in order to follow the dynamics further, a more complex 
Tij function might be numerically tested, holding the ef- 
fect of nonlinearities which are responsible for a stronger 
coupling strength. 

Concerning coherent phenomena from an experimen- 
tal point of view in extended systems, we should empha- 
size the dynamical features which are also shared with 
other systems in the route to weak turbulence: (i) the 
existence of critical domains with saturated widths (sta- 
tionary clusters in the ST/ZZ regime of size L c ), (ii) the 



existence of a cascade of bifurcations towards a regime 
(ST/ZZ) where the fundamental traveling modes split 
in both space and time, and (iii) the existence of local- 
ized patterns which are strongly nonlinear (ST/DW). We 
expect that this results may provide a new insight into 
spatiotemporal chaos and complex systems. 

Acknowledgments 

M.A. Miranda is grateful to P. Collet and R. Ribotta 
for fruitful comments. This work has been partly sup- 
ported by the Spanish Contract No. BFM2002-02011 
and No. FIS2008-01126, and by PIUNA (University of 
Navarra, Spain). M.A. Miranda acknowledges financial 
support from the "Asociacion de amigos de la Universi- 
dad de Navarra" . 



[1] V. K. Vanag, L. Yang, M. Dolnik, A. M. Zhabotinsky, 
and I. R. Epstein, Nature 406, 389 (2000). 

[2] I. Z. Kiss, Y. Zhai, and J. L. Hudson, Science 296, 1676 
(2002). 

[3] M. P. K. Jampa, A. R. Sonaware, P. M. Gade, and 

S. Sinha, Phys. Rev. E 75, 026215 (2007). 
[4] S. Rajesh, S. Sinha, and S. Sinha, Phys. Rev. E 75, 



011906 (2007). 

[5] G. Kozyreff, A. G. Vladimirov, and P. Mandel, Phys. 

Rev. Lett. 85, 3809 (2000). 
[6] K. Wiesenfeld, P. Colet, and S. H. Strogatz, Phys. Rev. 

Lett. 76, 404 (1996). 
[7] P. C. Hohenberg and B. I. Shraiman, Physica D 37, 109 

(1989). 



12 



S. H. He, H. B. Huang, X. Zhang, Z. X. Liu, D. S. Xu, [34 
and C. K. Shen, Phys. Rev. E 74, 057203 (2006). [35^ 
J. Gomez- Gardenyes, Y. Moreno, and A. Arenas, Phys. 
Rev. Lett. 98, 034101 (2007). [36 
A. C. Marti and C. Masoller, Phys. Rev. A 342, 344 
(2004). [37 
Y.Pomeau, Physica 23D, 3 (1986). [38 
A. T. Winfree, J. Theoret. Biol. 16, 15 (1967). [39 
A. T. Winfree, The geometry of biological time (Springer, 
1990). [40 
Y.Kuramoto, Chemical oscillations, waves and turbu- 
lence (Dover, 2003). [41 
D. Golomb, D. Hansel, B. Shraiman, and H. Sompolin- 
sky, Phys. Rev. A 45, 3516 (1992). ' [42 

J. Burguete, D. Maza, and H. Mancini, Physica D 174, 
56 (2003). [43 
J. Burguete, H. Mancini, and C. Perez-Gara'a, Europhys. 
Lett. 23, 401 (1993). ' [44 

D. Maza, J. Burguete, and H. Mancini, Int. J. Bifurcation [45 
and Chaos 4(5), 1353 (1994). 

R. Alvarez, M. van Hecke, and W. van Saarloos, Phys. [46 
Rev. E 56, R1306 (1997). 

E. Ringuet, C. Roze, and G. Gouesbet, Phys. Rev. E 47, [47 
1405 (1993). 

E. Ringuet, S. Meunier-Guttin-Cluzel, C. Roze, and [48 
G. Gouesbet, Int. Comm. Heat Mass Transfer 26, 477 
(1999). [49 
G. Gouesbet, C. Roze, and S. Meunier-Guttin-Cluzel, J. 
Non-Equilib. Thermodyn. 25, 337 (2000). [50 
L. Pastur, M. T. Westra, D. Snouck, W. van de Water, 

and M. van Hecke et al., Phys. Rev. E 67, 036305 (2003). [51 
L. Pastur, M. Henriot, and R. Ribotta, Phys. Rev. Lett. [52 
86, 228 (2001). 

A. Joets and R. Ribotta, J. Phys. 50, 171 (1989). [53 
L. Pan and J. R. deBruyn, Phys. Rev. E 49, 483 (1994). 

F. Giorgiutti, A. Bleton, L. Limat, and J. E. Wesfreid, [54 
Phys. Rev. Lett. 74, 538 (1995). [55^ 
A. J. Simon, J. Bechhoefer, and A. Libchaber, Phys. Rev. 
Lett. 61, 2574 (1988). [56 
L. Bellon, L. Fourtune, V. T. Minassian, and M. Rabaud, 
Phys. Rev. E 58, 565 (1998). [57 
P. Brunet and L. Limat, Phys. Rev. E 70, 046207 (2004). 
P. Brunet, J.-M. Flesselles, and L. Limat, Europhys. Lett. [58 
56(2), 221 (2001). 

C. Counillon, L. Daudet, T. Podgorski, and L. Limat, [59 
Phys. Rev. Lett. 80, 2117 (1998). 

P. Bot, O. Cadot, and I. Mutabazi, Phys. Rev. E 58, 
3089 (1998). 



C. Misbah and A. Valance, Phys.Rev. E 49, 166 (1994). 
R. J. Deissler and H. R. Brand, Phys. Rev. Lett. 72, 478 
(1994). 

R.J.Deissler and H.R.Brand, Phys. Rev. Lett. 81, 3856 
(1998). 

L.Gil, Europhys. Lett. 48, 156 (1999). 

P. Brunet, Phys. Rev. E 76, 017204 (2007). 

M. Miranda and J. Burguete, Phys. Rev. E 78, 046305 

(2008). 

F. Daviaud, J. Lega, P. Berge, P. Coullet, and M. Dubois, 
Physica D 55, 287 (1992). 

A. Kudrolli and J. P. Gollub, Phys. Rev. E 54, R1052 
(1996). 

S. Bottin, F. Daviaud, P. Manneville, and O. Dauchot, 
Europhys. Lett. 43, 171 (1998). 

D. Coles, In Mechanique de la Turbulence (CNRS Paris, 
1962). 

W. van Saarloos, Phys. Rep. 29, 386 (2003). 

R. Ribotta, A. Joets, and L. Lei, Phys. Rev. Lett. 56, 

1595 (1986). 

E. Plaut, A. Joets, and R. Ribotta, Journal de Physique 
III 7, 2459 (1997). 

M. Henriot, Ph.D. thesis, Universite Paris-sud, Orsay 
(2001). 

M. Dennin, D. S. Cannell, and G. Alhers, Phys. Rev. E 
57, 638 (1998). 

I. Oprea, I. Triandaf, G. Dalgelmayr, and I. Schawrtz, 
Chaos 17, 023101 (2007). 

M. Henriot, J. Burguete, and R. Ribotta, Phys. Rev. 
Lett. 91, 104501 (2003). 

Y. Tu and M. C. Cross, Phys. Rev. Lett. 69, 2515 (1992). 
S. Bottin, Ph.D. thesis, Universidad de Paris XI, Orsay 
(1998). 

B. Marts, K. Martinez, and A. L. Lin, Phys. Rev. E 70, 
056223 (2004). 

R.Ribotta, Phys. Rev. Lett. 42, 1212 (1979). 
K.Huang, Statistical mechanics (John Wiley & Sons, 
1987), 2nd ed. 

H. Daido and K. Nakanishi, Phys. Rev. E 75, 056206 
(2007). 

S. Mao, J. de Bruyn, Z. Daya, and S. Morris, Phys. Rev. E 
54, R1048 (1996). 

M. Miranda and J. Burguete, To be published in Int. J. 
Bifurcation and Chaos (2009). 

Y.Kuramoto, Proc. Int. Symp. on mathematical problems 
in theoretical physics, vol. 30 (Springer, New York, 1975). 



